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Abstract 

We demonstrate that the one-loop dilatation generator for the scalar sector of a 
certain perturbation of N = 4 Super Yang-Mills with fundamentals is the Hamiltonian 
of an integrable spin chain with open boundary conditions. The theory is a supersym- 
metric defect conformal field theory (dCFT) with the fundamentals in hypermultiplets 
confined to a codimension one defect. We obtain a /C-matrix satisfying a suitably 
generalized form of the boundary Yang-Baxter equation, study the Bethe ansatz equa- 
tions, and demonstrate how Dirichlet and Neumann boundary conditions arise in field 
theory, and match to existing results in the plane wave limit. 



1 Introduction 



M = 4 Super Yang-Mills theory in four dimensions has a number of fascinating properties. 
It is exactly superconformal. In the case of an SU(N) gauge group (to which we restrict 
ourselves here) it is believed to be precisely dual to Type IIB string theory in an AdS^ x S 5 
background with N units of five- form flux [T], j2], jHj, jlj. This duality becomes more 
tractable to study in the limit of excitations with large .R-charge, where the gravity dual 
geometry reduces to a maximally supersymmetric plane wave, as investigated by Berenstein, 
Maldacena, and Nastase In this regime it is particularly easy to see how a closed string 
arises as a collection of "bits," identified in the dual as fields in a gauge invariant single-trace 
operator. 

More recently it has emerged that an integrable structure is present on both sides of this 
duality. The planar one-loop matrix of anomalous dimensions for gauge invariant local op- 
erators has been shown to be the Hamiltonian of an integrable spin chain; the case of scalar 
operators was demonstrated by Minahan and Zarembo [E|, and their arguments were gener- 
alized to all gauge invariant local operators by Beisert and Staudacher [7j. Furthermore, on 
the string theory side the worldsheet Green-Schwarz theory was shown by Bena, Polchinski, 
and Roiban |8 j to possess an infinite number of nonlocal conserved quantities. The two 
manifestations of integrability were related by Dolan, Nappi, and Witten j^J; more recent 
work on integrability and AdS /CFT includes |l()j-|28j. 

One natural question to ask is whether the integrability of M = 4 SYM is present in other 
field theories. Systems closely related to M = 4 are the most likely to share in this special 
property; for other results on this question see [2H], (3D]- A particularly interesting feature 
to introduce is fundamental matter, associated with open strings in the gravity dual. 

One model with these properties is the defect conformal field theory (dCFT) formulated 
in |31j : for more work on defect field theories see [U2]-|IZ|- This system couples the four- 
dimensional fields and dynamics of TV = 4 SYM to M new three-dimensional hypermultiplets 
localized on a flat codimension one hypersurface, or "defect." These hypermultiplets trans- 
form in the fundamental representation of the three-dimensional part of the M = 4 gauge 
field, and break half the supersymmetry. There is only one free parameter, the bulk gauge 
coupling, and interestingly, the theory is exactly superconformal [31], |H21- For M <C N, the 
gravity dual is IIB string theory on AdS^ x S 5 , with the new ingredient of M D5-branes 
wrapped on AdS^ x S 2 [IE], [IS!- 1 Open string modes on the D5-brane are dual to gauge- 
invariant operators localized on the defect, which have the form of a number of bulk adjoint 
fields between a pair of fundamental defect fields; the defect fields thus play the role of the 
ends of the open string. 

In this paper, we consider the question of whether this dCFT is integrable, as M = 4 SYM 

1 For M ~ N the backreaction of the D5-branes must be taken into account, and the supergravity solution 
is unknown; but see [23!, EI]- 
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is. We study the sector of gauge- invariant, local defect operators composed of Lorentz singlet 
fields, generalizing the approach of [B] . After calculating the one-loop anomalous dimension 
matrix in the planar limit for this class of operators, we search for an integrable system 
with this operator as the Hamiltonian. This comes down to finding integrable boundary 
conditions for the SO (6) spin chain of [Oj, which can be codified in a matrix /C that solves the 
boundary Yang-Baxter equation, suitably generalized to incorporate boundary flavor degrees 
of freedom. 2 We find such a /C-matrix and demonstrate that the associated Hamiltonian 
indeed matches the one-loop dilatation generator of the dCFT, demonstrating integrability 
for our sector of operators. We conjecture that this integrability persists to the complete set 
of one-loop gauge invariant operators. 

Having obtained the integrable structure, one may use the Bethe ansatz to diagonalize 
the anomalous dimension matrix. Studying single excitations on the open chain and their 
reflections off the boundary, we demonstrate that they obey either Dirichlet or Neumann 
boundary conditions, matching the possible excitations of the dual open string. Furthermore, 
the periodicities and phase shifts conspire so that the the propagation of left- and right- 
movers on the chain is equivalent to that of only left-movers on an effective periodic closed 
chain of twice the length — the familiar so-called doubling trick of open string theory. The 
integrability of the system thus provides a window through which we may see more directly 
how open strings are realized within the dual dCFT. We also compare our results to the 
results of Lee and Park in the plane- wave limit [33], and find agreement. 

In section El we review the defect conformal field theory and calculate the planar one-loop 
matrix of anomalous dimensions. We solve the boundary Yang-Baxter equation to identify 
the integrable structure in section El and formulate the Bethe ansatz and consider single 
excitations in section HJ In section El we conclude. Field theory conventions are listed in an 
appendix. 

Note added. As this paper was being finalized there appeared work by Chen, Wang, and 
Wu [54J that obtains apparently similar results for a different deformation of M = 4 SYM 
with fundamentals. Interesting earlier work on open spin chains in gauge theories with 
matter appears in 55 -[60 . 

2 Field Theory Computations 

2.1 Review of Defect Conformal Field Theory 

The defect theory we study [2] is a descendant of Af = 4 Super- Yang-Mills theory with 
SU(N) gauge group. The M = 4 theory contains a gauge field A^, and adjoint Majorana 
spinors X a and real scalars X 1 in the 4 and 6 respectively of the SO (6) R-symmetry. The 

2 For work on classification of integrable boundary conditions without this generalization, see [S2], |53| . 
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dCFT introduces M additional 3D hypermultiplets propagating on the hypersurface x 3 = 0. 
The new fields break the total symmetry group from PSU(2, 2|4) down to OSp(4\A), where 
the 4D conformal group 50(4,2) is reduced to the 3D conformal group 5*0(3,2), the R- 
symmetry is broken SO(6)r — > S0(3)h x SO(3)y, and the supersymmetry is cut in half. 
The defect fields are M complex scalars q m transforming in the (2, 1) of 50(3)# x SO(3)y 
and M complex 2-component fermions \l/ a transforming in the (1, 2). Conventions and the 
action for TV = 4 SYM can be found in the appendix. 

The bulk fields are arranged with respect to the preserved 3D symmetry group into a 
vector multiplet {A k , P+A a , Xy, D^Xjj} and a hypermultiplet {A3, P_A a , Xjj, D$Xy }, where 
k = 0, 1, 2, A = 1, 2, 3, and I = 4, 5, 6. P+A and P_A with P± = (1 ± 7 5 7 3 )/2 indicate the 
splitting of a 4D Majorana spinor into two 3D Majorana spinors. The Xh and Xy scalars 
transform as (3, 1) and (1, 3) of SO{3)h x SO(3)y, while the A a are arranged into a (2, 2) 
field X arn . 

The defect fields couple directly only to the bulk vector multiplet. The total action for 
the theory is that of TV = 4 SYM together with 

S3 = Skin ~\~ Sy U k + Spot j (1) 

S kin = j,Jd 3 x((D k q m yD k q m -z^ a p k D k ^ a ) } (2) 

Syuk = jtJd 3 X (l^ a P + \ am q m - Kr^maP^ + * a <T&X$* a ) , (3) 

S P ot = j,Jd 3 x {jTX$X$r + uL IJK iT°LX 3 H X%f) (4) 

+ ^ / d 3 x(q^aUD 3 Xj I )q n + \5(0)TT( ( ra I mn q n ) 2 ) , 

with P^g = d k q — iA k q and similarly for The total theory has only one coupling g, which 
is exactly marginal 32J. We have written (fT|) - for M — 1, which is sufficient for our 
purposes; for M > 1, dual to multiple D5-branes, flavor indices appear on q and \1/ in an 
obvious way. 



2.2 Anomalous Dimensions for Scalar Operators 

We are interested in local, gauge-invariant defect operators constructed out of Lorentz scalar 
fields that are "single-trace" in the sense that one line of color runs through the entire 
operator. These have the form 

= ^ m j u ... jL , n q m X j K..X^q n , (5) 

with the obvious sum over color indices implied. We can think of each operator as an open 
spin chain consisting of vector spaces W& x V\ x . . . x Vl x Wp, with Vi in the 6 of 5*0(6), 
(3 in the 2 of SO(3) H , and a in the 2. 3 

3 The 2 is equivalent to the 2, but we write it this way to avoid worrying about epsilon symbols. 
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We are interested in the one-loop planar contributions to the matrix of anomalous dimen- 
sions for these operators. We calculate this by evaluating the correlation function 

T = (q n ,(z p )X^(z L ) . ..X^(z l )q m '(z a ) O(0)) . (6) 

Only nearest neighbor fields along the chain will interact to this order, reducing the contri- 
butions to effective correlators of two fields with a two-site piece of the operator O. 

First consider the interaction between two bulk scalars X ta and X %b . Interactions that exist 
purely in the M = 4 theory were calculated for the case of the closed spin chain jB] (see also 
, [SI], |M| and include the nontrivial interactions of gauge boson exchange and the 

4-point X interaction, as well as contributions of the one-loop wavefunction renormalization 
of each X field. In principle, additional interactions involving the defect fields can contribute 
as well. However, one can quickly convince oneself by looking at the action JIJ that there 
are no nontrivial defect interactions between 4 X's at lowest order. Furthermore, although 
diagrams with defect fields contributing to (XX) do exist, they are forbidden on symmetry 
grounds from leading to a divergent piece in the X self-energy |31j . 

Consequently, the contribution of two X fields to the matrix of anomalous dimension is 
unchanged from the pure M = 4 case. As a result, all the interactions in T not involving a 
defect field match the contributions in the case of the closed chain. (This also implies that 
the anomalous dimension matrix for closed chains of X's in the defect theory takes exactly 
the same form as in M = 4 SYM.) The defect interactions are localized to the ends of the 
chain and will have the interpretation of integrable boundary conditions added to an already 
integrable system. 

2.3 Differential Regular izat ion and Bulk Interactions 

Although they are well known [HJ, we include the calculation of the XX interactions to 
establish our conventions. We perform our computations in position space using differen- 
tial regularization jMj. The presence of the defect favors a position space regulator, and 
differential regularization is particularly convenient for extracting anomalous dimensions. 

The essential philosophy of differential regularization is to replace amplitudes containing 
position-space distributions too singular to have a Fourier transform by regularized distribu- 
tions differing only at the singular point; this replacement is equivalent to adding local coun- 
terterms. The regularized amplitudes have a well-defined Fourier transform but depend on 
an arbitrary mass scale; differentiating with respect to this scale produces the usual Callan- 
Symanzik equations for the regulated amplitudes. As we are interested only in anomalous 
dimensions, we keep only logarithmic divergences, and discard quadratic divergences. 

Consider first the amputated X self-energy T x = (X l (z\)X^ (z 2 )) ■ There are three one- 
loop diagrams — fermion loop, gauge boson loop, and scalar tadpole; the last is a purely 
quadratic divergence and hence irrelevant for us. The fermion and gauge boson diagrams 
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give 

= g 2 N5^(4 - 2)d»A 12 d»A 12 = g 2 N5 ij -^^- , (7) 

l 47T ) Z 12 

where we discarded a quadratic divergence going like A 12 D A i2 ~ b(z\ 2 )jz\ 2 . Regularizing 
according to the formulae in the Appendix, we find 

M ^7 T x = ^^n5(z 12 ). (8) 



DM 87T 



Since /3(g) — for our dCFT, the Callan-Symanzik equation is particularly simple: 

DM 



M^r 7 + 2 lx (g))T x = 0. (9) 



7^) = Vr- (10) 



Using the amputated zeroth-order result = — (1/2)5 AB n 5(z 12 ), we find 

g*N 

8vr 2 

To find the contribution of a pair of X fields to the anomalous dimension of the defect 
operator, we compute the correlator of two X fields with the appropriate part of O(0): 

r 4 = (x h (zi)x i2 (z 2 ) x^x»(o)) . (ii) 

Amputating the external legs associated with X 11 and X 12 , regularizing, and keeping only 
the divergent parts, we arrive at 

m Jm v ^ = 9 ^ 6{zi)6{z2) ( (7 ^- + (73 ^- + (Ts) ^) ' (12) 

where the tensor structures coming from the gauge boson exchange, 4-point X interaction, 
and self-energy between X n and X 71 are 

(^)& = <M, = 2%% (T E )^ = -4#$ . (13) 

Only one of the two possible self-energies is included, so that the other may be included in 
the interaction of the next pair in the chain. Correspondingly, we include only one factor r y x 
from the Callan-Symanzik equation for the total correlator T: 

( M ^ + 7x + 7o ) r4 = °- (14) 

Using this equation, the zero-loop expression Tf^ = N5 i J \5 l ^ l 5(zi)5(z 2 ), and the result (fTUj) 
for 7x, we obtain 

Ma = {2m - + shi2 Kn} ■ as) 

Summed over all nearest neighbor pairs, ()15j) forms the complete one-loop planar dilatation 
operator for a closed chain [Hj. 
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2.4 Defect Interactions 

We turn now to the defect interactions. One necessary ingredient is the q self-energy. Defining 
T q = (q m (yi)q n (y2)) , we find 4 diagrams at one-loop order: a \I//A fermion loop, a d^Xn loop, 
a gauge boson loop, and a q tadpole; the last diagram does not contribute to the logarithmic 
divergence. The contributions to the amputated self-energy are 

r£i /A = 2g 2 N5™Tr(S 12 h2), (16) 

= -/iV5™ (2d k V 12 d k A 12 + 2A l2 V 2 V l2 + \V l2 V 2 A l2 ) , (18) 
leading to the total logarithmic divergence: 

after regularizing according to the formulae in the appendix. The amputated zeroth-order 
propagator is — 5™V 2 8{i)i 2 ), and hence the Callan-Symanzik equation implies 

(^ + 27^))r. = - 7, = 0. (20) 

Next we consider the interaction at one end of the chain. We calculate 

T d ^(X\z 1 )q m '(y 2 )q m X^0)). (21) 

A universal, flavor-blind interaction is gauge boson exchange, which leads to 

g 2 N 2 S^S ij (A 12 d k A 1 d k V 2 - \A,V 2 V 2 A l2 + ^A^V^An - ±V 2 d k A 1 d k A 12 ) . (22) 

If one pulls out derivatives to obtain divergence-free total derivatives plus a term with all 
derivatives on V 2 which contains the divergence, you get 

r« ^ = total derivatives + ^^^Wfo)- - M— r« p = |^C><^i)<% 2 ) .( 23 ) 

There is another diagram with a scalar interaction qqXX, which is different depending on 
whether the two X scalars are X v or X H . For X % = S l A Xy, X J = 5 3 B Xy, the q and X scalars 
are charged under different 5*0(3)8 and the flavor interaction must be trivial. We find 

rgW = -^sz's^siz^v, - ^^ v = -^'^WW- (24) 

In the case X 1 = 5}Xjj,X^ = 8 3 K Xjj the flavor interaction is nontrivial: 

rgL = -^fs(z 12 )(ieuKrtJV 2 A 2 -> Jlflr^ = -^(z^WMl) .(25) 
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The Callan-Symanzik equation for the complete correlator T is 

+ L lx + 2 lq + 70) T = , (26) 

where V and 70 are understood to be matrix-valued, while the other terms are proportional 
to the identity. The matrix 70 can be separated into the parts acting in the bulk of the chain 
on two X's, and the boundary pieces. There are L — 1 pairs of neighboring X's, each of which 
contributes to the matrix as (fTo^) . Notice, however, that this does not take into account one 
factor of 7x in (jZEj) , as well as one contribution of the X self-energy to T; this is because 
the open chain has no X lL X n interaction. Because these contributions are proportional to 
the identity in flavor-space, we can insert them anywhere; we find the convenient choice to 
be to include half the contribution of each with each boundary term. 

The boundary interaction at the left end of the chain is then determined by 

+ f+7, + »)r, + iMA r (.) = c. (27) 

We obtain 

{lo)it: = (2C'$ + 2ie IlJlK a^ m ) , ( 70 )™;i: = (utfs£) , (28) 
and analogously for the other end of the chain, 

(7o):::K = - 2^^,) , { l0 y±i = , (29) 

where delta-functions in the unwritten indices are understood. 
In summary, we find the matrix of anomalous dimensions to be 

lo = 16^2 I ^ haa+1 + {2hl + 2S ^ + (2/l/3 + 2Sl ^ ) ' (30) 
haa+l — K aa +i + 2J aa+ i — 2P aa +l , (31) 
where subscripts on operators indicate the vector spaces they act on, and 




JSJaJb 

iaib 


= S jajb 5- ■ 
u u i a ib 1 


pjajb 
iaib 


U tb la ' 


TjaJb _ 

iaib 


U la U lb ' 




qIttl 
°Jn — 


K 

—ie IJK a nm , 


oAm 
°Bn 


xms:A 
— °n°B 1 


Qlm 
°Bn — 


cAm 
°Jn — 


0, 


Qfhl 
D nJ 




QfhA 


cm, zA 

= ° n 6 B) 


Qfhl 
°rLB _ 


QfhA 
^hj — 


0. 



(32) 



A check of our computation is that the anomalous dimensions should vanish for protected 
operators. In [31], the chiral primaries of the dCFT were determined to be 

r°<LX%X I i...X I H ^\\ (33) 

where parentheses denote total symmetrization and tracelessness among all the triplets of 
SO (3). One may indeed verify that these operators are annihilated by 70 as given in (j3~Uj) ; 
in fact each h a ^ a+ i, as well as 2/ s i + 2Sai and 2I L/3 + 2Sl/3, separately gives zero. 
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3 The Boundary Yang-Baxter Equation 



Having calculated the one-loop planar dilatation generator for the class of defect operators 
(jSJ), we want to learn whether it is the Hamiltonian of an integrable open spin chain. 

It was shown by Minahan and Zarembo in [0] that the operators composed of closed 
chains of X fields possess this kind of integrable structure. This was done by identifying an 
7£-matrix, in this case the S'0(6)-invariant 

K l2 {u) = X - \u{u - 2)/ a2 - (u - 2)P 12 + uK 12 ) , (34) 

where 1 and 2 label the two vector spaces acted on by the operator, and I, P, and K are as 
defined in (j32J). These satisfy the Yang-Baxter equation 

n 12 (u)n 13 (u + v)n 23 (v) = n 23 (v)n 13 (u + v)n l2 {u) . (35) 

If we define the transfer matrix as the trace of the monodromy matrix, 

t(u) = Ti a T a (u) = Tr a K al (u)K a2 ■ ■ -7l aL {u) , (36) 

where a is an auxiliary space and the dependence of t{u) and T a {u) on the vector spaces 
1,2, ... L comprising the chain is understood, the Yang-Baxter equation guarantees that the 
transfer matrices commute for arbitrary arguments: 

[t(u),t(v)] = 0. (37) 

The expansion of the transfer matrix in powers of u then generates an infinite number of 
conserved quantities, starting with the momentum and Hamiltonian. The identification of 
the one-loop dilatation operator for closed chains as the Hamiltonian associated with (|3*4*j) 
demonstrated the integrability for that system. 

The analogous method for introducing integrable boundary conditions was formulated in 
[65 ; for a pedagogical summary see jHE] or section 3.5 of \§7\. In addition to an integrable 71- 
matrix determining the dynamics for the bulk of the chain, one introduces operators JCf(u) 
acting on either end of the spin chain. These satisfy the boundary Yang-Baxter equations, 
or BYBs, 

n X2 {u-v)Kl{u)n l2 {u + v)K^{v) = K^{v)n X2 {u + v)K^{u)n l2 {u-v) (38) 
K 12 (v - u^t^i^TZ^i-u -v- 2ij)!Ct t2 (v) = Kt t2 {v)K 12 {-u -v- 2i 7 )/C^ 1 {u)K 12 {v - u) . 

Here, ti refers to taking the transpose on the ith vector space, and 7 is a parameter charac- 
teristic of a given 7?.-matrix such that 

K%(u)7l\\{-u - 217) = \{u) (39) 
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for some scalar function X(u) (for our case it is easy to see that 7 = 2i). We can then form 
the open chain transfer matrices 



t(u) = Ti a K, + a {u)T a {u)lC-{u)T^{-u) . 



(40) 



If (JHSJ) and (|3*8)l are satisfied, these commute, 



[t{u),t(v)) = 0, 



(41) 



and thus generate the family of commuting charges characteristic of an integrable system. 

We must generalize this formalism slightly, from matrices /Cj~ and /Cj that act only on 
the last "ordinary" vector space in the chain, to matrices /C^ and /C7g that also act on the 
degrees of freedom W & and Wp living on the boundary In this case the resultant t(u) acts 
on the vector space W a x V\ x • • • x V L x Wp. 

The structure of integrability is largely unaffected by this modification; the extra indices 
"go along for the ride" in the BYB, and (j41j) still holds. Given the quantum numbers of Wp 
and Vl-i the matrix K~(u) must take the form 



where rh, n are indices in Wg, I, J are indices in the 3 of SO(3)h, and A, B are indices in the 
3 of SO(3)y- Furthermore, if we have a /C - that satisfies the first BYB equation (jHEJ), then 
(/C/ - )' 1 (m) = K-i (2 — u) satisfies the second BYB equation; however, this process will give 
matrices acting on a spin chain whose boundaries are both in the anti-fundamental of SU{2). 
We can easily translate one of the boundary actions to act on the associated fundamental, 
and we will do so later. 

Thus what we need to do is find what functions f(u), g(u), and h(u), if any, satisfy the 
first BYB. Writing out the BYB with indices, we want 



Ufliu - v)(K-)^(u)U k e X(u + v){K-)f^v) = (/C-)™>)7^> + v)(JC-)Z(u)US:(u - „).(45) 



There are six independent indices in this equation: 21,^2,^1,^2 and m,p. The first four can 
each take a value in either SO(3)h or ,SO(3)y, allowing for a total of 16 equation types that 
need to be satisfied. 

It is not hard to see that both sides will vanish whenever three of the indices ^1,22,^1,^2 
belong to one of the 5*0 (3) 's and the fourth to the other 5*0(3). This leaves us with 8 
equations to check. Out of concern for the reader, we will give only the equations and their 
solution, omitting details. The case where all indices belong to the SO(3)v turns out to 
be satisfied independently of the form of the f(u), g(u), and h(u), as do the cases with 




h{u) 5™5i , 
(IC-)Z A (u) = 0, 



.K 
ran i 



(42) 
(43) 
(44) 



Mi e SO(3) H or i 2 ,£ 2 e SO(3) H . 
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The same two conditions arise from either £±,£2 G SO(3) H or ii,i 2 G SO(3) H : 

(u + v) [(u - vf + 2] [h(v)f(u) - h(u)f(v)} + (u-v) [{u + v) 2 + 2] [h(u)h(v) - f{u)f{v)\ 

= 2{u + v + l)(u - v){u + v - 2)g{u)g{v) , (46) 

{u 2 — v 2 ){u — v — 2)[h{y)g{u) + h{u)g(y)} + (u + v){u — v — 2)[h{v)g{u) — h{u)g{v)] 
-(u 2 -v 2 )(u + v- 2)[g(u)f(v) - g(v)f(u)} -(u-v)(u + v- 2)[g(u)f(v) + g(v)f(u)} (47) 
= -i(u + v + l)(u - v)(u + v - 2)g(u)g(v) - 3{u 2 - v 2 )g(v)(h(u) + f(u)) . 

The cases «i, £2, G SO{3)h and 22,^1 G SO(3)h give two more equations: 

(v-u)[f(u)f(v)-h(u)h(v)-2g(u)g(v)\ = (v + u)[h(u)f(v) - h(v)f(u)\ , (48) 
(v - u) [f(u)g(v) + f{v)g{u) - ig{u)g(v)} = (v + u) [h(u)g(v) - h(v)g(u)} . (49) 

The case with all indices in S0(3)h, which is the most labor-intensive, gives two final 
equations: 

ig(u)g(v)(u-v)[2u + 2v-l] = 2vf(v)g(u)-2uf(u)g(v), (50) 

ig{u)g{v){u — v){u + v — 2)[2(u + v){u — v — 2) + u + v + 1] = 
2vg{u)f{y)[u 2 - v 2 + 2} - 2g(v)f(u)[u(u 2 - v 2 ) - 3{u 2 - v 2 ) + 2u}- (51) 

3g(v)[f(u) + h(u)](u 2 -v 2 ) . 

We find that 

f(u) = 2u 2 -u + 1, g(u) = 2iu, h{u) = -2u 2 + u + 1 (52) 

satisfy these conditions, up to an overall function of u that we fix by requiring /C~^(0) = Ip L . 
(This overall ambiguity corresponds to reshuffling the basis of the infinite set of mutually 
commuting operators inside i{u).) We therefore obtain the result 

(£")#(«) = (2u 2 -u + + 2uie IJK a« fi , (53) 

= (-2u 2 + u + m5^ (/C-)^(«) = (/C-)J[^(«) = 0, 

and this implies 

(IC + )™J(u) = (2u 2 -7u + 7)5^5 I J + (2iu-U)e IJK a^ fi , (54) 

(lC + )™£(u) = (-2u 2 + 7u-5)5^, (/C + )™iH = (/C + )™/(«)=0. 

Remembering that we want to deal with a spin chain having one boundary in the fundamen- 
tal, we need to determine, given the above Kr acting on the anti-fundamental, a Kr that 
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acts on the fundamental. For SU(2), if ip is in the fundamental, then a 2 ip is in the anti- 
fundament 
we obtain 



fundamental; therefore, we should have 1C = o~ 2 /C a 2 . Recalling that a 2 a K a 2 = —(a K ) T , 



K 



= (-2« 2 + u + 1)6*6$ , (^-)S(«) = (^-)^(«) = 0. 
Now, we want to form the transfer matrix 



(55) 



tin) = Tr n 



K.i a {u)Tl al {u) ■ ■■K aL (u)iC~Ju){K La {-u))- 1 ■ ■ ■ (KU-u))- 1 



Since (as is straightforward to verify) 



K 



-U 



-7Z ab (u) , 



we find 

i(u) = Tin 



{A-u 2 ) L {l-u 2 ) L 



(A-u 2 )(l-u 2 ) 
£ta(u)n al (u) ■ ■■K aL (u)jC- f3 (u)K aL (u) ■ ■■TZ al (u) 



(56) 



(57) 



(5? 



We are interested in expanding i(u) = Y, n inU n , and identifying the Hamiltonian as t\. 
Making use of the relations 



n ab (0) = P ab , j-K ab {<S) = l - [-21 ab - P ab + K ab ] 



(59) 



and the definitions in (}32|) . we determine that the lowest term in the expansion of i(u) is 

"4 L 



t = t(Q) = Tr a 



^I^Sa(0)-Pal • • • PaLhpPaL ' ' ' P a l 



Tin 



/C; 



61. 



(60) 



The next term is slightly more complicated. Breaking it into pieces, we find that 



and 



Tr a [dJC+a (0)K al (0) ■ ■■K aL (0)jC- /3 (0)K aL (0) ■ ■ -K al {0) 

Tr a \dJCta (0)P a l • • • PadapPaL • • • Pall = Tr a \d u 1C+ a (0) 



Tr 

■ LL o 

Tr 

i T r 

2 a 



£+ (0)d u TZ al (0) ■ ■■TZ aL (0))C- p (0)TZ aL (0) ■ ■■TZ al (0) 
/C+(0)^ al (0) • • • 7^(0)^(0)^(0) • • -d u n al (0) 

— 2I a l — P a \ + K a i)P a2 ■ ■ ■ PaLlafiPaL ' ' ' PqX 



+ 



K+a (0)P al • " • PaLhaPaL " " " P^i^hl - P al + 



al i 



Tr„ 



/C+ (0)(-2/ al - P al + K al ) = \2S al - 13/, 



(61) 



(62) 
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The terms with the derivative acting on lZ a £ for £ £ 2, L all behave similarly: 



Tr a [£+ (0)ft al (0) ■ ■■d u K al {Q) ■ • -^(0)^(0)^(0) ■ • -72 al (0) 

Tr a [/C+ a (0)72 al (0) • • • 72^(0)^(0)72^(0) • •• ^72^(0) • • -72^(0)' 
1 



2 
1. 



:Tr r , 



^(O)Pal 



-2U - P ai + K a 



Tr a JC^ a (0)P a i ■ ■ ■ PaLlfiaPaL 



-21, 



• PaL^apPaL " ' " Pal 
Pat + K a t)---P ai 



+ 



+ 



Tin 



/C+ a (0)P al • ■ ■ P^_! [-2P a , - J a , + A^] P a ,_! ■ • ■ P 



ol 



Tr„ 



/C+(0) (-2P £ _ W -J/_ w + K, 



e-u 



Q(-2P U - I. u + K t 



Finally, we have a term of the form 



Tr a [/C+ a (0)72 al (0) • • • 72 aL (0)«9^(0)72 ai (0) • • -72 ai 



Tr„ 



K.%* (0)P„i • • • Pa L d u K-J0)P aL • • • P 



ol 



Tr„ 



Tr„ 



£+ (0)1 ^^(0) = 6 



KiMPai---d u Kl^) 



We combine these expressions to give 



(63) 



(64) 



Pal 



h = -13/ai + 12^i + Sti (-12P^+i - 6/ wl + 6K i4+1 ) - 6I/3L + 12S L(3 
and adding a multiple of the identity to the result, we find 

h + (25 + 18L)/ = 6 \£^h u+ x + (2/ s i + 25 a i) + (2I pL + 2S Lp ) 



(65) 



(66) 



which is indeed proportional to the matrix of anomalous dimension 7e> found in (|H0jl. demon- 
strating that this operator is in fact a Hamiltonian for an integrable system. 



4 The Bethe Ansatz 

Having seen that the planar one-loop matrix of anomalous dimensions (fHU|) is an integrable 
Hamiltonian for the spin chain, we may use the techniques of the Bethe ansatz to find the 
eigenvectors and eigenvalues of this matrix. 

To do so one begins with a Bethe reference state, corresponding to a certain operator, and 
systematically creates excitations above this ground state. These excitations correspond to 
impurities in the chain. The "Bethe ansatz" is to assume that the eigenvectors consist of 
waves of impurities (or "spin waves") propagating along the chain with certain momenta 
hi] integrability assures that scattering of these waves factorizes into products of two-body 
scattering where the magnitude of each momentum is separately conserved. 
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The Bethe ansatz requires for consistency one equation for each excitation. For the closed 
chain, these equations are of the form 



e ikiL = llS ji (k j ,k i ), (67) 

where the S- matrix Sji(kj,ki) is the phase shift acquired when the i th excitation passes 
through the j th excitation; the equations (|6*Tj) can be thought of as the requirement that 
the phase of an excitation as it travels all the way around the chain and comes back to 
itself is unity. In the special case of non-interaction, Sji = 1 and (|H7j) reduces to the usual 
quantization of momentum on a circle. 

In the case of systems with Lie group symmetry, extra excitations are introduced that 
correspond to changing the orientation in group space of the spin wave. For example, consider 
the general formula given in jHj, eqn. (4.35): 

/ u q>i + ia q ■ w/2 \ L _ Ugj - u q<j + ia q ■ a q /2 tt tt u g,j ~ u g',j + ^g • Qg'/ 2 ^ 

\Uq }i -ldq-w/2) j*. U q>i - U q J - %OLq ■ OL q /2 ££ q ^ U q>i - Uq>j ~ lOL q • Ct q '/2 ' 

Here the u q ^ are parameters characterizing excitations, taking the place of the kf, i labels 
the excitation as before, while q reflects the fact that the excitation can be associated to any 
of the simple roots a q of the algebra, w is the highest weight vector of the representation of 
the group that lives at each site. For a fundamental at each site, we will have w = w l , the 
first fundamental weight, which has the inner product with simple roots a q ■ w 1 = 8 q . 
Thus for an excitation u q i with q — 1, the left hand side is 



»q,i 

L 



ui,j + i/2 
«i i - i/2 



(69) 



which allows us to match the equation to (JS7|) by using the relation between ui ti and the 
quasi- momentum ki (jEj, eqn. (4.39)): 

= -i log 777: • (70) 

Uij — i/z 

Note that the momenta and energy depend only on the u\ t i excitations. For the u q> i, q 7^ 1 
excitations, the left hand side of (p)Hj) becomes trivial, but the right hand side constrains the 
possible u q j. One should think of the u\ t i as equivalent to the quasi-momenta k^ and creating 
an actual excitation with energy and momentum, while the other u q ^ do not change the 
energy or momentum but instead change the Lie group quantum numbers of the excitation. 

The open chain Bethe ansatz equations generalize (|H7jl in an intuitive way. In the open 
case, instead of making one complete circuit of the closed chain, picking up phases Sji for 
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each interaction, a given excitation passes one way across the chain, reflects off one boundary, 
passes back the other way across the chain, and reflects off the other boundary: 

e 2 iki L = B ^ k .) B ^ k .) -Q ( Sj .( kjj ki)Sji(kj, -hi)) . (71) 

Here B\ and £>2 are phases generated by reflecting off the boundary, while the Sji are the same 
S-matrix terms as in the closed chain; the interaction between the spin waves is ignorant of 
the boundary conditions. We get each Sji twice because we pass both ways along the chain, 
with the relative momentum reversed between the two interactions. 

So we see that if one already knows about the closed chain, the Sji are known, and 
determining the open chain Bethe ansatz equations comes down to finding B\ and 82- In 
our case the Sji are determined by (J55J) for group 50(6), as was studied in [§]. The reflection 
coefficients B± and £>2 involve only the excitation in question and the boundary; whether there 
are other excitations is irrelevant. They can hence be determined by considering the case 
with just a single excitation. We do this below. 

4.1 Bethe Reference State 

We use as our Bethe reference state a particular chiral primary of the dCFT. Defining 

Z = X\ + iXjj, we take 

\0) L = q 1 ZZ...Zq 2 , (72) 

which results from (|33)l where for each a = 1 . . . L + l one has chosen I a = +, in the standard 
basis where J3 is the Cartan generator of 50(3)//. This state has jn = L + 1, jv = 0. 

4.2 Single Impurity Sector — SO(S)v Fluctuation 

We now consider the case where any single Z field in the chain is replaced by another X 1 
such that jn = L. A There are two ways we can do this: we can excite one of the 50(3)y 
scalars, or we can excite Xfj. 

First consider an SO(3)v excitation; since the Bethe state is invariant under 50(3)y, any 
of the three scalars are equivalent, and we will let W be any of Xy, X v , or X v . We define 
the state |W(a;)) of the spin chain as follows: 

\W(x)) = q l ZZ...W...Zq 2 , (73) 

Note that replacing a Z by Z actually leads to jn = L — 1; this mixes with two X 3 impurities and hence 
will not be included in this section. 
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where W has replaced Z at the x site. We then make a standard Bethe ansatz for the 
energy eigenstates \W(k))\ 

\W(k)) = £ f(x)\W(x)) , f(x) = A(k)e ikx + A(k)e- ikx . (74) 

x=l 

Unlike the closed chain ansatz, where we need consider waves propagating only in one direc- 
tion, this case requires, because of reflection, that both directions be present. Integrability 
underlies the assumption that the magnitude k of the momentum will not change, either in 
interacting with another excitation or in reflecting off the boundary. 

Now we wish to determine the action of the Hamiltonian on the position eigenstates 
\W(x)); for convenience we work with H defined as 70 = (g 2 N / 'I6ir 2 ) H. Consider first the 
case where 1^ 1, L; the boundary parts 27^1 + Rm, 21 + 2Sl/3 of H act on the same fields 
q\Z and Zqi as in the chiral primary, and hence give zero. The same is true for hjj+i acting 
on a pair ZZ. Hence the only nonzero contributions to the Hamiltonian are 

h x , x+1 \W(x)) = (K XtX+1 + 2I XtX+1 - 2P x , x+1 \W(x)) = + 2\W(x)) - 2\W(x + 1)) , (75) 
h x -^\W(x)) = (K x . liX + 2I x . ltX - 2P x ^, x \W(x)) = + 2\W(x)) - 2\W(x - 1)) , (76) 

leading to 

H\W{x)) = 4\W(x)} - 2\W(x - 1)) - 2\W(x + 1)) . (77) 

Now consider the states |W(1)) and \W(L)). Nonzero contributions are 

h lt2 \W(l)) = 2\W(1)) - 2\W(2)) , (2/ al + 2R m )\W(l)) = 4\W(1)) , (78) 
h L -i,L\W(l)) = 2\W(L)) - 2\W(L - 1)) , (2I Lf3 + 2R Lp )\W{L)) = A\W{L)) . 

Thus we have 

H\W(1)) = 6\W(1))-2\W(2)), H\W{L)) = 6|W(L)> - 2\W(L - 1)) . (79) 
Now we demand that the states \W(k)) be energy eigenstates: 

(H - E(k))\W(k)) = 0, (80) 

implying the relations 

= (A-E(k))f(x)-2f(x + l)-2f(x-l), 2<x<L-l, (81) 
= (6-£(fc))/(l)-2/(2), = (6-E(k))f(L)-2f(L-l). 

Using the first relation we can deduce 

E(k) =4(1 -cos A;). (82) 



15 



We may then solve for A(k) and A(k) as follows. All the relations (|81|) are consistent with 
the single result (J%2|) only if the exceptional second and third equations in (|%T|) take the same 
form as the first, with auxiliary quantities /(0), f(L + 1) introduced. This follows only if 



/(0) = -/(l), /(L+1) = -/(L). 

If we plug this into the Bethe ansatz (JUJ), these equations imply 

A(£;)(l + e* fc ) = -A(k){l - e~ ik ) , 
A(fc)e ifcL (l + e ifc ) = -i(A;)e- tfcL (l+e- ifc ), 

which can be seen to have a solution only if 

e 2ikL _ j 



(83) 



(84) 
(85) 



*6) 



This is the Bethe ansatz equation for this impurity. We see that in fact the boundary 
interactions are trivial and an elementary quantization of the quasi-momentum k is implied: 



B X = B 2 



k 



rcn 



17) 



Meanwhile one finds the solution A(k)/A(k) = — (1 + e lk )/(l + e tk ). The energy eigenstates 
()74j) can hence be written in terms of the integer n: 



\W(n)) 



a 



L 

n Sil1 

x=l 



nn 



x - 1/2) 



\W(x)) 



up to some normalization a n , with the associated anomalous dimension 

9 2 N 



lW{n) 



Alt 2 



COS 



nn 



19) 



We see that the boundary conditions imposed on the spin waves are that they go to zero at 
the "end points" a half-step past the last links on the chain; these are Dirichlet boundary 
conditions. 

This is precisely what we expect from the dual string theory point of view. In the string 
dual, this operator corresponds to an open string living on a D5-brane; and oscillations in 
the three SO(3)y directions normal to the brane must obey Dirichlet boundary conditions 
at the endpoints. Moreover, the periodicity (|86|) with no B factors tells us that the left and 
right movers on the chain can be recast as just left moving excitations on a closed chain 
of length 2L; this is nothing but the familiar "doubling trick" of open string theory. We 
now turn to a discussion of the X 3 excitation, which should be associated with Neumann 
boundary conditions. 
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4.3 Single Impurity Sector — X 3 Fluctuation 

We now consider an X 3 impurity, which we will abbreviate as just X. This lowers the jn 
charge to jn = L without turning on a jv charge. As before we define states with a single 
impurity 

\X{x))=q 1 ZZ...X...Zq 2 , (90) 

with X replacing Z at site x. However, these states are not a closed set under action of the 
Hamiltonian. We may also create states with ju = L, jy = by flipping the spin of one of 
the q fields on the end of the chain, but leaving all bulk excitations as Z: 

|<?i) = q x ZZ ...Zq-i, \q 2 ) = q 2 ZZ . . . Zq 2 . (91) 

We assume the energy eigenstates to be of the form 

\X(k)) = ]T g(x)\X(x)) + a\ Ql ) + (3\q 2 ) , g(k) = B(k)e lkx + B(k)e~ ikx . (92) 

x=l 

The analysis then follows as before. Requiring \X(k)) to be an energy eigenstate 

{H - E{k))\X{k)) = (93) 

leads to the equations 



= 


(4 


- E)0 + 2g(l), 




(94) 


= 


(4 


-E)g(l)+A(3-2g(2), 




(95) 


= 


(4 


-E)g{x)-2g{x-l)- 


2g(x + l), 2<x<L-l, 


(96) 


= 


(4 


-E)g(L)-2g(L-l)- 


■ 4a , 


(97) 


= 


(4 


-E)a-2g(L). 




(98) 



The energy has the same form as in the previous subsection, as can be deduced from the 
generic (2 < x < L — 1) equation 

E(k) = 4(1 — cos A;) . (99) 

Again we proceed by attempting to cast the other equations in this same form. For the 
second and fourth equations, this can be done using the definitions 

g(0) = -2(3, g(L+l) = 2a. (100) 

The first and last equations also share this form if we introduce the fictitious 

g(-l) = g(l), g(L + 2)=g(L). (101) 
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These relations imply for the ansatz ([92)1 



B{k){e- ik - e ik ) 
B(k)e ikL (e 2ik - 1) 



B{k){e~ ik - e ik 
B(k)e 



ikL Q ^—2ik\ 



The first obviously forces B(k) = B(k). The second is then consistent only if 



3 2ifc(L+l) 



k 



nn 



L + l 



(102) 
(103) 



(104) 



Once again the Bethe ansatz equation indicates that the boundary reflection coefficients 
are trivial: B\ = B% = 1. The chain is effectively one link longer here, because of the 
participation of the q fields in the spin wave. The energy eigenstates can be written 



L+l 

\X(n)) = b n Y^ cos 

x=0 



irnx 
L + l 



\X(x)) 



where we have defined 



|X(0)> = -(1/2)|&>, |X(L+l)> = (l/2)| ft > 



and the anomalous dimension is 



lW(n) 



4vr 2 



cos 



nn 



L + l 



(105) 



(106) 



(107) 



We notice immediately that for these excitations, the boundary conditions are Neumann. 
Since X corresponds to a fluctuation in a direction in which the brane is extended, this is 
again exactly what we expect from the dual string point of view; the possibility of flipping 
the spin on the q fields plays a fundamental role in allowing the motion of the endpoints. 
Again the left and right movers together behave as a single closed chain, this time of length 
2(L + 1). 

We conclude that the integrability of the system allows us to see with particular clarity 
the emergence of open strings from the defect local operators. Oscillations of the string are 
nothing but spin waves propagating on the chain, reflecting with boundary conditions that 
are equivalent to trivial transmission onto a "doubled" closed chain. 



4.4 Plane Wave Limit 

The plane wave limit of the gravity dual to the dCFT was studied by Lee and Park [33J. 
They considered a Penrose limit where the defining null geodesic wound around the S 2 C S 5 
on which the D5-brane is wrapped; the associated light-cone vacuum \0,p + ) is just our 
Bethe reference state, the chiral primary \0)l, with (fip + a') 2 = J 2 /(g 2 N), J = L. They 
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furthermore studied single- impurity excitations, and it is interesting to match these to our 
results. 

Lee and Park identified an open string excitation in a Dirichlet direction having oscillator 
number n with the operator 

1 J v^sinf— ) 

4' D |0,P + ) - ^jE N J/2 + i ^Z l X 3 Z J - l q 2 , (108) 

while a Neumann direction excitation took the form 5 

1 J \/2cosf— 1 

a]l N \0,P + ) - NJ/ lr ^Z l X*Z J ~ l q2 . (109) 

They identified the anomalous dimension for both cases as 

_ irg s Nn 2 

Using the relation Aitg s = g 2 , we find that the large- L limit of our results indeed agrees with 
the conclusions (jl08j) - (jll0|) of for the plane wave limit. 



5 Conclusions and Open Questions 

The presence of integrability inside four dimensional quantum field theories is unexpected 
and fascinating. We have seen that the sector of open chain operators constructed from 
scalar fields in the superconformal defect deformation of Af = 4 Super Yang-Mills possesses 
an integrable structure, taking the form of integrable boundary conditions added to the 
structure of M = 4 SYM alone. It is natural to hypothesize that this integrability holds for 
all operators in the dCFT. Given that theories closely related to M = 4 are also integrable, 
one may speculate on how broad a class of theories possesses an integrable structure, and to 
what extent integrability is a useful concept for understanding other, non-conformal gauge 
theories, for example those that confine. 

The existence of the open chains also implies that in the gravity dual, open string boundary 
conditions on D-branes can be integrable. It is natural to suspect that an OSp(A\A) Yangian 
will survive in the open Green-Schwarz worldsheet theory. It would be interesting to verify 
the presence of this symmetry, and to look for other D-brane configurations that preserve 
integrability. 

Finally, we have seen in the interactions of the spin waves in the Bethe ansatz how prop- 
erties of open strings like Dirichlet or Neumann boundary conditions and the doubling trick 

5 The authors of mention that operators analogous to 1)91(1 with q spins flipped could mix with (|1()9(1 . 
but do not include them explicitly. 
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arise in the field theory dual as an open string is assembled from "bits" of quantum fields on 
the spin chain. In this sense integrability, like supersymmetry, provides a crutch to help us 
see the workings of duality It would be interesting to use the Bethe ansatz to diagonalize 
more complicated states so as to learn more about mesonic-type operators and their proper- 
ties under duality, and in greater generality to learn more about how strings are assembled 
from field theory constituents. 
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Appendix: Field Theory Conventions 

Our index conventions are as follows, fi, v and k, I are 4D and 3D Lorentz indices. We denote 
the 6 of SO(6)r with it decomposes under SO(3)h x SO(3)v as 6 — > (3,1) © (1,3), 
labeled with I, J and A, B, respectively. The 4 of SO(6)r is and it decomposes as 
4 — > (2, 2); doublets of SO(3)h are denoted m, n and doublets of SO(3)v are a, b. 

We work in Euclidean space rotated from mostly-minus signature, where the 7-matrix 
conventions are those of [31], and are such that Majorana spinors are real. 

We write our adjoint fields in matrix notation, and color indices are implicit throughout. 
The action for M = 4 SYM is 

S A = j,J d 4 xTr {\F^ - i\ a YD^\ a + D^X a D»X a - \{X\X b ] + C a ap \ a [X\ A' 3 ]) (111) 

with the covariant derivative D^X a = d^X a — i[A^, X a ] and similarly for A a . 

The Yukawa tensor C°p in the 4 x 4 -> 6 can be written explicitly as follows. Split 
the fermions as \ A = A a — 1,2,3, and A = A 4 , and assemble complex scalars as 
$ A = (X$ + iXfi)V2 with Xy 2 ' 3 = A 1 - 2 ' 3 and X 1 /' 3 = A 4 ' 5 ' 6 . Then 

J d 4 x Tr C^X a [X a , A^] = (112) 

i2V2 J d 4 xTr (\[Lx A , $ A ] - X A [R\ ® A \ - ±e AB cX A m C + R$ C ),x B ]), 
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with L = (1 + 7 5 )/2, i? = (1 - 7 5 )/2. 

The A Q are arranged into X am transforming in the (2, 2) of SO(3)h x S , 0(3)y: 

Aim = A5 im - ix"V^ • (113) 

We use z^ to indicate 4D coordinates and decompose them into 3D coordinates y k and 
transverse x = x$. The massless 4D and 3D scalar propagators in Feynman gauge are 

(X( Zl ) X(z 2 )) = ^A 12 = ^ 2 j~v 2 > (AXzi) A^ 2 )> = g V A i2 , 

(A(z x ) A> 2 )> = -S^ = V^A 12 , (g(2/i) g(ife)> = 2>u s , (114) 

where we have suppressed the color delta functions; the factors of 1/2 arise for the 4D fields 
because of the matrix notation we are using. Since A is Majorana, we also have non-vanishing 
(A(zi) A (22)) and (X(zi) X(z 2 )), determined by A = A T 7°. 
We write the 3D fermion \1/ in a 4D notation, as a four-component spinor satisfying 

P + ^ = #, P_$ = 0, P ± = |(l±7 5 7 3 ) , (115) 

where 7 5 = —ijjj 2 ^ is the chirality matrix. The propagator for \1/ is then 

(*( Zl )*(z 2 )) = ( 1 k P + ) d k V 12 = S X 2 ■ (116) 

The projection matrices P± obey 

[^±,71 = 0, ^P ± = P^ X , Pl = Pl = l, P+P- = 0. (117) 

Propagators for X im can be written in terms of (|114j) as 

(X im (z 1 )\ nj (z 2 )) = 2{X(z 1 )X(z 2 ))S ij 5 mn , 

(X im (zi) X jn (z 2 )) = 2(X(z 1 ) X(z 2 ))e ij e mn , (118) 
(X mi (zx) X nj (z 2 )) = 2(X(z 1 ) X(z 2 )) 

In differential regularization we use the replacements (neglecting quadratic divergences): 



1 

? 


1 log z 2 M 2 
* 1° z 2 ' 


(119) 


1 

1? 


1 log z 2 M 2 
* 32°° , 2 ' 


(120) 


1 

1 2/ 1 3 


v ^og M\y\ 

\y\ 


(121) 


1 

M 5 


l^ 2 hgM\y\ 

6 \y\ 


(122) 



where □ and V 2 are the 4D and 3D Laplacians, respectively. 
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